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AGISSILAOS ATHANASSOULIS AND THIERRY PAUL 

Abstract. We present several results concerning the semiclassical limit 
of the time dependent Schrodinger equation with potentials whose regu- 
larity doesn't guarantee the uniqueness of the underlying classical flow. 
Different topologies for the limit are considered and the situation where 
r£ two bicharacteristics can be obtained out of the same initial point is 

emphasized. 
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1. Introduction 

The state of a quantum system is described by a positive, trace class 
operator D £ , which evolves in time under the Heisenberg-von Neumann 
equation 

is§- t D%t)= [-^A + V,D%t) 
D £ (t = 0) = £>§. 

Both the positivity and the trace of Z? £ are preserved in time |20j ; observ- 
ables, also described by operators, can be measured against the state. For 
example, the value of observable G 6 , associated with operator G £ ', at time t 

l 
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is given by the trace Q £ {t) = tr(G £ D £ (t) ). The term density matrix is often 
used for a positive, trace-class operator. 

The Heisenberg-von Neumann equation is closely related to the Schrodinger 
equation. By virtue of the singular value decomposition, the density ma- 
trix can be decomposed to a sum of pure states, i.e. rank-1 orthonormal 
projectors, 

D£ = yZ X m\ u m) ( ^ n \ 



m 



(the A^'s are time- independent in the standard context of an external, time- 
independent potential, i.e. if V does not depend on D £ , or t. Here we limit 
our scope to such problems). Then, each of the if^'s satisfies a problem of 
the form 

i£ -m u m = ( ~T A + V ( x )) u£ m, 



n 2) "^dt^rn — y 2 

<(< = o)a 2 

The small parameter e is called Planck's constant, and when it tends 
to zero one expects the system to behave in a way approaching classical 
mechanics; i.e. motion under Newton's law in a force field with potential 
V. This is a singular asymptotic limit, and several techniques have been 
developed for its study; one is based on the Wigner transform (WT) [TTJ Q21 
V. 



The WT W £ {x, k, t) is defined as 
(1.3) W £ {x,k,t) = W £ [D £ (t)](x,k) = J e- 27Tiyk K £ (x + e^,x-e^,t)dy 



where for each t, K £ (t) is the integral kernel of the operator D £ (t). The WT 
is also closely related to the Weyl symbol of D £ , and the Weyl calculus of 
pseudodifferential operators in general, see e.g. |lUj for more details in that 
direction. It evolves in time under the well-known Wigner equation, 

d t W £ (x, k) + 2nk • d x W £ (x, k) + 



(1.4) +IRe ife 27TiSx V(S)W £ (x,k- f ,t)dS 



0. 



W £ (t = 0) = w { 



o - 



an equivalent reformulation of (1.1). In this context, it can be shown that 
observables can be measured by G £ {t) = J G £ v (x,k)W £ (x,k)dxdk, where 
G £ v (x,k) is the Weyl symbol of the operator G £ . (However this in general 
is not a Lebesgue integral, but only makes sense as a Cauchy PV one). 
In the case of a pure state, D £ = \u £ )(u £ \, the WT can be written as 

(1.5) W £ (x,k,t) = J e- 2niyk u £ (x + e^,t)u £ (x-e^,t)dy. 



SEMICLASSICAL LIMIT FOR ROUGH HAMILTONIANS 



3 



As a quadratic transform on an oscillating function u £ , of small typical scale 
of oscillation e, it has been widely used in signal processing, see e.g. [9] and 
a multitude of references therein for more on that direction. 

Under appropriate conditions the WT allows for a very natural and com- 
pact description of the semiclassical limit. Indeed, it has a physically mean- 
ingful limit as e tends to zero, while in general the operator D £ itself does 
not. The limit (in an appropriate sense; see below for more precise statement 
and reference), is called the Wigner measure, 

W £ (t) -» W°(t), 

and is a probability measure that evolves in time under the Liouville equation 
of classical statistical mechanics 

d t W° + 2irk ■ d x w° - ±d x V ■ d k w° = 0, 

(1.6) 

W n ° = UmWZ. 

u e^O u 

Observables are now described by bounded continuous functions, which can 
be measured against the classical state, G°(t) = f G°(x, k)dW°(t). The 
formal probabilistic interpretation e.g. for the expected position of a particle 
X(t) = J xdW°(t) should be intuitively clear. 

Thus the correspondence between classical and quantum physics for the 
same system is highlighted in a very intuitive way, in terms of models that 
are, at least formally, similar and directly comparable. Convergence issues 
are not always painless or completely understood, as we will see in more 
detail. There is also the well-known caveat that W £ (t) is in general neither 
a measure of bounded total variation, nor non-negative e.g. as a distribution. 
Finally asking that the limit problem itself should be well-posed isn't always 
trivial. 

Note at this point that considering initial data with 

W £ {0) 5(x-x ,k-k ) 

is particularly natural: they stand after all for classical particles with posi- 
tion xq and velocity ko. 

It was shown in [14J (Theoreme IV. 1) that, under appropriate conditions 
on the initial data (which allow for practically all physically meaningful 
choices), as long as the potential is smooth enough, V G C , the WT con- 
verges (in weak-* sense, and up to extraction of a subsequence) to a weak 



solution of the limit problem (1.6). It was immediately commented that, for 
V £ C 1 \ C ' , this does not yield an effective computation of the semiclassi- 
cal limit in general, but only a partial description: even after the extraction 
of a subsequence (so that there is a unique lira Wn n ) , in general there are 

n 

several weak solutions for the Liouville equation that W £n (t) could converge 
towards (because the Liouville equation with C 1 \ C ' potential and mea- 
sure valued initial data has existence, but not uniqueness of weak solutions. 
This can be seen e.g. using the Peano Theorem). Therefore some kind of 
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selection principle would be needed to supplement the result, if we were to 
be able to identify the Wigner measure corresponding to our subsequence 
in e of initial data. This, in a nutshell, is the context for the present work. 
Our goal is to understand better the mechanisms that may lead to loss of 
uniqueness, and develop tools which, at least in some cases, can resolve 
them. 

We are particularly concerned with initial data that concentrate VF £ (0) — 1 
5(x— xq, k—ko), and our findings here apply, under appropriate assumptions, 
to such data. The main assumption (other than, essentially, C 1,e potential 
regularity), is, roughly speaking, a non-concentration condition on the flow 
associated with the Liouville equation. 

Before the statement of the present results, some reference to related 
works is warranted. The problem of extending semiclassical asymptotics 
to non-regular potentials (e.g. worse than C 1 ' 1 ) has been studied by many 
authors - and for good reasons: some of the simplest, most natural systems 
(such as a gravitational or electrostatic force field) are described by poten- 
tials that are not C , 1, or C . In most cases the conclusions make use of 
a particular form of singularity (e.g. Coulomb/piecewise smooth potential 
etc) as opposed to a general smoothness class [UJ [13], and/or need some 
additional, non-trivial condition (non-concentration, non-interference etc) 
[171 E] ■ Another type of results is for whole random populations of initial 



data [2] (which in particular can be even weaker than a result applying to 
"almost all" initial data; i.e. there might be no way to just choose "almost 
any" initial datum at t = and keep track of it; the conclusions apply the 
population as a whole). Finally, it must be mentioned that in the study of 
the semiclassical limit in nonlinear problems (e.g. the Schrodinger-Poisson 
equation, a very natural problem), a lot remains to be understood. One 
of the main difficulties, comes precisely from the lack of regularity of the 
potential (which now depends on D £ (t)), and the subsequent lack of well- 
posedness for the corresponding classical problem for measure- valued initial 
data; see e.g. [IS EES [2]]. 

More specifically, semiclassical limits with rough potentials were consid- 
ered recently in [2\ where the hypotheses on the potential are, roughly speak- 
ing, that it has a generic part with BVi oc gradient, plus possibly a repulsive 
Coulomb part. The setting involves a random population of initial data, and 
it is shown that the corresponding population of solutions at a later time 
tends weakly to the push-forward by the Liouville equation of a population 
of measures. The classical propagation problem for the random family of 
measures is shown to be well-posed by virtue of the Ambrosio-Di Perna-Lions 
theory [H [7] . In [S] similar core ideas were used with deterministic mixed 
states. In that context, an averaging condition forbids initial data concen- 
trating e.g. to a delta function in phase-space, a fact consistent with the 
weak limit strategy and the fact that the flow is not defined everywhere. The 
present work therefore concerns less general potentials than those treated 
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in [21 18] , (still giving rise to ill-posed classical dynamics of course) , but can 
deal with initial data concentrating to point-supported measures. 

The possible ill-posedness of the limit equation can originate not only from 
lack of smoothness of the potential, which we focused on for our discussion 
so far, but e.g. from working on a torus instead of Euclidean space [To] . 

Finally, let us present a concrete example. Already in fT2] , families of 
initial data were constructed so that their Wigner measures at later times 
oscillate between the several solutions of the Liouville equation. The one- 
dimensional case study used for that, was essentially equivalent to the fol- 
lowing: 

Consider the potential 

(1.7) V(x) = -\x\ 1+e ■ (3(x) 

with G (0,1), C°° 3 /3 = 1 on and (3{x) = for \x\ > 2, and 

initial data such that Wq — 1 S(x — 0, k — 0). Then one can easily check that 
there exist, for t small enough, two families of bicharateristics starting near 
the origin at t = 0, namely 

(1.8) (X±(t),P±(t)) = (±c t»,±^- 1 ), 

with v = and Co = ^ - ^ ^ 

In this case it is possible to find semiclassical limits where one wavepacket 
(a (5-function in classical phase-space) splits in two wavepackets. (In this 
particular example it happens immediately, i.e. Vt > 0, or not at all; different 
examples can be constructed where it happens after an e-independent time 
i*). The mass of each wavepacket depends on information which can be 
easily extracted from the quantum problem, but is lost if we take the limit 
as usual in a straightforward way. 



The first of our results deals explicitly with the case study from [T3] men- 
tioned earlier, to provide a clear illustration to the selection principle at 
work. Our second Theorem gives a general result concerning the approxi- 



mation, in strong topology, of the Wigner function of the solution of (1.4) 
with the solution of a "smoothed" Liouville equation. Our third results goes 
back to the weak approximation, but with explicit remainder estimates. 

2. Main results 
We will use a space of test-functions introduced in |14j : let 

(2.1) A = {4> e C (R 2n ) | fsu P \^ K [^(x,k)}\dK<^}, 

equipped with the norm 

(2.2) IHU = j sup\T k ^ K [<j)(x,k)]\dK. 



(i 



A. ATHANASSOULIS AND T. PAUL 



Our first result will focus on the case study presented in the introduction; 
a family of initial data which concentrates at the origin at time t = and 
splits into two separated wave packets at time t > 0. This is only an example 
among several others presented in section [5] and the more general results of 
the paper are given by the Theorems |2.2| and |2.3| below. 



Theorem 2.1. Let V given by (1.7) and 



7+39 , 1+S 1-9 , / 2\ n o x 2 +k 2 

(2.3) W£{x,k) = \^r W {\^Tx,\^k) * I - J e^^— 

with A = log(l), w G H 2 nL°°nL 1 ) w(x,k) ^ 0, suppw C {|x| 2 + |/c| 2 < 1}, 
J w(x, k)dxdk = 1. 

Then 3T > s-uc/t f/m* /or a// f G [0,T], f/ie solution W £ {t) of (|L4j) 
converges in weak-* sense in A! to 

(2.4) W°(t) = c + ,5(X + (t),P + (t)) + c_<5(X-(t),p-(t)) ; 



with (X ± (t),P ± (t)) given by (|1.8l) and 



c± = / w(x,k)dxdk. 

±x>0 



- 



Remark: The fact that the initial data contain a smoothing by (f ) e 
ensures that Wq corresponds to a density matrix, i.e. a positive trace-class 



operator (see e.g. Exemple III. 7 of [H], or Lemma 6.7) 



In Remarque IV. 3 of [13], one can find a construction of a family of nor- 
malized wave functions uf concentrating at the origin, and whose Wigner 
functions, after extraction of a subsequence, converge to a Dirac mass cen- 
tered on (X + (t), P + (t)). The same construction is obviously possible in 
order to get a family u~ whose Wigner functions, after extraction of a sub- 
sequence, converge to a Dirac mass centered on (X~(t), P~(t)). Defining 
Dq = ~j^(\uf)(uf\ + \u~)(uj\) one can check that the Wigner function of 
Dq will follow the conclusion of Theorem 2.1 The present construction is 



different, and is not obtained through compactness arguments: there is no 
need of extraction of subsequence, and the scaling property is explicit (it is 
given by an implicit diagonal argument in a two scale sequence in |14j ) . The 
construction of [H] is based on the more-or-less explicit understanding of the 



flow around the singular point, while our approximation (i.e. Theorems 2.2 



and 2.3 below) is built around a non-concentration condition, which is rela- 
tively easy to check for other problems, whether the flow can be computed 
explicitly or not. (This may be increasingly important in high dimensions). 

Let us also point out the similarity with the long-time behaviour around 
a regular separatrix, studied extensively in [18]. In the regular case, the 
quantum wavepacket leaves the fixed point in "infinite" (i.e. asymptotically 
long) time, as does the classical flow; here the classical flow leaves the fixed 
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point immediately, and this behaviour (under additional assumptions, of 
course) is preserved by the quantum system. 

We turn now to the more general results of this paper. 

In the sequel we will consider potentials V that satisfy the following as- 
sumption: 



Assumption 1. 



S 



2 



(2.5) / \V(S)\ — ^ dS <oo, 

and moreover there are constants C > 0, 9 G (0, 1) such that for m G {0, 1, 2} 

Vl^aO^+oo: / \V(S)\\S\ m dS ^ ^^(b™- 1 - 9 -a™- 1 - ) 
\S\e(a,b) 

This is closely modeled after V(x) = C\x\ 1+9 ; indeed it is easy to check 
that the aforementioned potential satisfies this condition. In section [5j we 
will also see some other relevant types of singularities (also generated from 
| x | 1+8 in some sense). 



Denote 

(2.6) V(x) = (^\ A J e-^ x - x '^V(x')dx'. 

Theorem 2.2. Let us suppose that Assumption [7] holds, and there exist 
T > 0, 5 G (0, §^0), 7 > and approximate initial data [Wq] G H 2 (M. 2n ) 
such that 

(2-7) \\m}-W£\\ L 2= (\\W£\\ L2 ) 
and the solution of 

(2.8) d t p + 2irk'd x p-—d x V-d k p = 0, 
with initial condition p(t = 0) = [Wq] satisfies 

(2.9) \\p(t)\\ H i = 0(e- 5 \\W£\\ L 2) 
uniformly on [0, T] . 



Then the Wigner function W £ (t) of the solution of (1.4) satisfies, uni- 
formly on [0, T], 

(2.10) \\W £ (t)-pf(t)\\ L i=0(e K \\WS\\ L 2 + ||W ^-[^]|| i2 ) = o(||^|| i2 ). 



where p\ is the solution of (2.8) with initial datum p\{t = 0) = Wq ant 

2+6 



k = min{ 7^4^ — 1, 2TD ~ & }• 
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Remark: The non-concentration condition. It is helpful to discuss 
the motivation behind this type of assumption. Ideally, one would like the 
solution of the problem 

dtp + 2nk ■ d x p - ±-d x V -d k p = 0, 



2tt' 

to have H 2 regularity in the sense ||/?(i)||#2 ^ e ct \ \ Wq \ \ H 2 for some C that 
is either e-independent, or grows extremely slowly in e. Of course, when 
working with low-regularity potentials, this is in general not possible. So 
one looks for a weaker version of regularity, which would still be sufficient 
here; it turns out that the property 

3[W$] G H 2 {R 2n ), \\\W§] - Wfi\\ L 2 = o{\\W^\\ L 2) such that, if 

d t p + 2wk ■ 8 x p - ±8 X V-d k p = 0, p(t = 0) = [WJ] 
then ||/9(i)||#2 is not too large 

is an acceptable surrogate. 

We call this a non-concentration condition, because of one of the simplest 
ways to check it. Assume that there is a "small" set S C I 2 ", such that 

the derivatives of V up to order 3 outside S are "not too large" in e. Now 

set B = |J cfi £ _t(S), where cftf is the flow associated with the regularized 

te[o,T] 



pit = 0) = W$ 



Liouville equation (2.8). B contains all the "dangerous" part of phase- 
space; away from B the original regularity requirement holds automatically. 
A "concentrating flow" could pull too big a part of phase-space through 
the original "small bad neighbourhood" S in time t G [0,T]. On the other 
hand, a "repulsive singularity" , in general tends to push the flow away, thus 
keeping the pre-image of S small enough. The approximate initial data [Wq] 
can be then constructed by restricting the initial data on B . 

Remark: Checking the assumptions for concrete problems. Con- 
crete examples for which the assumptions (in particular the non-concentration 
condition) can be checked to hold, are examined in section[5] The techniques 
used can be easily extended to similar problems. It must be noted that the 



condition for Theorem 2.2 is seen to be satisfied for a generic selection of 
initial data, as long as the rate of concentration satisfies certain constraints 
(i.e. the profile that concentrates can be any function with a given regu- 



larity; see e.g. Lemmata 5.1| 5.2). In contrast, for the stronger Theorem 



2.3 (see next Remark for the difference), in general we need to restrict the 
shape of the wavepacket with additional assumptions for the same problems; 
Lemmata 15.31 15.41 



2.2 



is an 



Remark: L versus measures, and convergence. Theorem 
1? asymptotic approximation result: two functions, of constant L norm in 
time (also equal to one another and large in e), are relatively close in L 2 . 
It is not, however, a convergence result, since in general neither function 
has a limit in L 2 . Ideally, we would like to conclude that, since they're 
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(relatively) close, both functions converge (in some appropriate sense) to the 
same limit, typically a (sum of) delta function(s). However, the (relative) 
I? approximation is not by itself sufficient to ensure that if e.g. pf(t) — 1 
5(x — x(t), k — k(t)), then also W £ (t) goes to the same limit. Is there some 
additional information in the problem that we could use to come to that 
conclusion? Or is there a physical issue here, of 1? norm really "missing" 
the solution? For example, is it really possible that in a problem with 



splitting (as in Theorem 2.1) the conclusion of Theorem 2.2 holds, in the 
limit p\ "goes to the right", but at the same time W £ (or half its mass) 
"goes to the left"? This is a question that we cannot really answer at the 
moment. One possible way to go is by strengthening the non-concentration 



condition, and that is pursued in Theorem 2.3 below. This yields full control 



over the convergence (in an appropriate sense, see also Corollary 2.4), but 
at the cost of excluding many interesting, natural problems. 



Theorem 2.3. Assume that the hypotheses of Theorem 2.2 hold and p\{t) is 
as in Theorem 2.2. Moreover, in addition to the condition (2.1) we assume 
that 



(2.11) 



WW 



'oWl 2 



o(l), 



'o\\l 2 



o(e K ) (k being the same as in Theorem 2.2), and the initial 



data Wq corresponds to a density matrix, i.e. it is the Wigner function of 
a positive trace- class operator Dq, triD^) = 1 Ve > 0. 

For <j) E A denote 

E £ {4>) = (w £ {t)- P \{t),4>). 

and, for brevity, by F(e) the bound of equation (2.1(fy , namely 
F(e) = C(e K \\W^\\ L 2 + \\[W£] - W^\\ L .) = o(l), 
for some large enough constant C > 0. 



(2.12) 



Then, V</> G A, t ^ 0, the following estimate holds: 

n 

|^(^)| ^ F(e)M\ A + -^ T) VW)+ 



where x 6 C°°(M 2n , [0,1]) is any cutoff function, xi z ) = 1 ^\ z \ < 1> x( z ) 
V|z| > 2, and xm(z) := x(Mz). 



Remark: I? versus measures, and convergence. If a generic "nice" 
function concentrates, f £ — 5(x,k), and — [/ £ ]||li = o(l), then it follows 
that [f 6 ] — 5(x,k) as well. In other words \\\W( 



general holds in the examples where Theorem 2.2 does) would better suited 



Wq||^i = o(l) (which in 



to control concentration. We are forced however to use the much stronger 



L 2 = o(l) due to technical reasons, and in particular the fact 
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that the Wigner equation doesn't respect the L 1 norm in general (or the 
total variation of measures, or even the A' norm away fro m th e positive 
cone which also scale like L l ). The assumptions of Theorem 2.3 should be 



considered pessimistic, and probably can be relaxed. In any case, the precise 



convergence result is summarized in Corollary 2.4 below 



Remark: Size of the error. The last term in rhs of equation (2.12) can 
be further estimated if there is more information on cfi, and in particular its 
decay in the k variable. 



Corollary 2.4. If the assumptions of Theorem 2.3 are satisfied and ||[Wg]||£i 
is uniformly bounded in e, then 



UmE £ ((j)) = 
e^O 



e A, t ^ 0, 



and therefore the Wigner function W £ (t) has a weak-* limit in A' if and 
only if p\ has a weak-* limit, and the two are equal, i.e. 

(2.13) 3W°(t) : lim(W £ (t),<p) = (W°(t),<p) V0 6 A, t € [0,7], 

e^O 



if and only if 
(2.14) 



lim(p\(t), <p) = (W°(t), 4>) V<f>eA, te [0, T]. 

£^0 



In other words, under the assumptions of Theorem 2.3 one has to solve 



equation (2.8) in order to compute the propagation in time, p\ (t), and then 
take the limit e -> 0. That is, unlike what happens with regular potentials, 
the semiclassical limit and the propagation in time by (a regularization of) 
the Liouville flow no longer commute. 



p?(o)- 

£^0 

¥ 

W°(0) 



~p\{t) 

£^0 

W°(t) 



Figure 1. In some cases, there is not enough informa tion in 
M /0 (0) to determine W°(t), i.e. the Cauchy problem (1.6) is 
ill-posed. 



Organization of the paper: Proofs of the main Theorems are in section 
|4j Checking the assumptions of the Theorems for various concrete problems 
is done in section [5j Certain auxiliary results we use are collected in section 

El 
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3. Definitions and Notations 
The Fourier transform is defined as 
(3.1) f(k)=T x ^ k [f(x)]= J e~ 2 ™ kx f(x)dx. 

For compactness, we will use the following notations: 

i J e 2niSx V (S)W(x, k - f)dS 



TYW = \R& 

c s 



2F _1 



JV(S)W(X - S, K) 



sin(ireSK) 



dS 



}-d x V -d x W 



2tt 



2-kF 



-l 

X,K^x,k 



f V(S)W(X - S, K)S ■ KdS 



The Sobolev norms of order m on phase-space will be defined as follows: 

ll/l|w m .P(K 2 «) = ^ \\dxdkf\\LP(R 2 ")i 
|o|+|6|<m 

where of course a and b are multi-indices of length n each. Moreover, 
H m (R 2n ) = W m ' 2 {R 2n ). 

4. Proofs 

4.1. Proof of Theorem 12. 11 In this section we will show that Theorem 



2.1 follows from Theorems |2.2[ |2.3| The latter are proved in the next sub- 
sections. 

Before getting into more detail, let us comment that A as used in the 



statement of Theorem 2.1 is an effective upper bound. It is easily checked 



that if we replace A by A' ^ A, y = o(l), the Theorem still holds. 
Remark on notation: In this section we use the notation R = A - 3. 

Now denote 

(4.1) fo{x,k) = \^Tw(\~^x, ATI - *;), 



where of course A = log(^), i.e. Wq = /q * (f )™ e 2?r " 
that 



. One should note 



(4.2) 



\W? - m\ L2 < V^\\m\m = 0(e-*-*>) Vr? E (0, i) 



(The approximation follows from a standard observation on smoothing op- 
erators, see e. g. |4 | for a proof. Moreover, since by construction the constant 
k of Theorem 2.2 is « < |, the error \\W® — /q||l 2 can indeed be dropped 
without loss of generality, when invoking Theorem 2.3 - and Theorem 2.2 
through it). 
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In Lemma 5.3 it is shown essentially that Theorem 2.3 applies if n = 1, 



V(x) 



in {|x| < 1} (with a smooth cutoff outside of that, the 

k 



h t ). where S \ 



5.3 



specifics of which are irrelevant) and W§ = 5 x 1 S k 1 w( 

R~2~ ; ^ = R~ . (We avoid duplicating parts of the proof of Lemma 
here. Moreover, it is completely straightforward to see that the assumptions 
not checked explicitly there also hold). Therefore Corollary 2.4 applies, i.e. 
it suffices to find the limit in e of p\ . 



Claim: With [Wq] defined as in equation (5.15), we have 



(4.3) l|Wo £ -[Wo £ ]|Ui= (l). 

Proof of the claim: We will work separately for the two terms 



m - [WoIHli < \m - /oIIli + 11/0 - ra^- 



First of all, it is obvious that 



\ft 



W^dxdk 



\W^\dxdk 



0(e c 



|(£C,fe)|>2 l(a:,fc)|>2 

Moreover, using the Jensen inequality, one sees that 



\(x,k)\<2 



W§\dxdk < 



/ dxdk f |/ £ 

\(x,k)\<2 \J\(x,k)\<2 



W$\ 2 dxdk ^ 



and therefore 



iw? 



/pll/ 



o(l). 



Moreover, 



(4.4) 



w{ 



\x\<CR 



5 X ' s k 



) | dxdk 



J \w{x,k)\dxdk^C5- 1 R = o(l). 



P2, 



The proof of the claim is complete. 

This means that p\ and p\ (the evolution under equation (2. 
Wq respectively) are interchangeable for our purposes, since 

{p\A) = {pIA) + {p\ 

and 

\{p\-pIA)\^\\px- 

(We used the obvious bound 
eq. d2lb ). 



of [W$], 



l°° < \\ws-[w$]\\ L i mu . 

Loo ^ ||</>||.a; recall that A was defined in 
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In particular, this means we can always work with |x| > CR, since p\ stays 
outside of {|ac| < CR} by construction. This allows certain ODEs we will use 
to be well-posed. Without loss of generality we will work for x > 0, k > 
(working with k < makes no difference other than the opposite sign in 
the explicit expression for K(0) in equation (4.5), the case x < follows by 
symmetry). Denote by X(t),K(t) the unique solution of 

X(t)=2rK(t), K(t) = ^(X{t))\ 

(4-5) 

X(0) = CR, K(0) = \X(Q)\ 1 ¥, 
i.e. the only branch of the level set of {2n 2 k 2 + V{x)} = in {x > R A k > 

The interesting property of this trajectory is that, unlike what happens 
with regular potentials, it leaves zero in finite time (more generally, i.e. 
without restricting to {x > R A k > 0}, it reaches and leaves zero in finite 
time) . 

Claim: Any characteristic of the Liouville equation (2.8) starting in {x > 
R A k > 0} Pi suppfWjf] converges to X(t),K(t) defined in equation \4- 
Then it readily follows that 

Xx>o(x, fc)[W £ ] ->> J Xx>o(x, k)w(x, k)dxdk 8{x - X(t), k - K(t)). 

Repeating the argument for x < gives a limit of two delta functions leaving 
zero along different trajectories. 

Proof of the claim: Denote by Xi{t),K\{t) the solution of 

X x {t) = 27r#i(i), K^t) = -±d x ViX^t)), 

(4.6) 

X 1 (Q) = x , Ki(0) = fc 0; 

where of course (xo,ko) £ {x > R A A;>0}n supp[Wq] are as in the 

statement of the claim). Then equation (4.6) can be recast as 

(4-7) 

Xi(t) = 2-KK^t), Kx{t) = -^V^Xrit)) - ±(d x V{X 1 (t)) - V^X^t))), 

X 1 {0) = x , K 1 (0) = k . 

Here we use the estimate \d x V(Xi(t)) - V x (Xi(t))\ < CR^e'i = o(l) (see 
Lemma 6.2 with V'(x) in the place of /). It is clear now that any such 
characteristic converges to 

X 2 (t) = 2vrK 2 (t), K 2 (t) = -^V x (X 2 (t)), 

(4.8) 

X 2 (0) = x , K 2 (0) = k . 
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To conclude observe that trajectory Xz{t), K%(t) is squeezed between two 
level sets of 27t 2 k 2 +V (x) converging to each other (namely between {2n 2 k 2 + 
V(x)} = V{CR) and {2-n 2 k 2 + V(x)} = ViR 1 ^)). 



Essentially the same analysis applies to the example of Lemma 5.4 as well 



4.2. Proof of Theorem 2.2[ It will be helpful to recall the main objects we 
are going to use here. The WT for this problem, W = W £ (x,k,t), satisfies 
the well-known Wigner equation (1.4) with initial data Wq. Moreover, W( 
and p\ are defined as follows 

( d t Wf + 2irk ■ d x Wf + T/Wf = 0, 



(4.9) 
and p\ 
(4.10) 



W{ (t = 0) = W§, 



d t p' 2 + 2irkd x p e 2 +T Q v p £ 2 = 0, 
p*(t = 0) = [W§], 

where [Wq] for now is simply assumed to exist (and have the properties 



found in the statement of Theorem 2.2 . Such approximate initial data will be 
constructed, for concrete examples, in Section[5]). Finally, recall that p\ was 
defined as the solution of equation (2.8) with initial data p\e{t = 0) = Wq. 



We partition the proof as follows: 

\\W £ - pl\\ L 2 < \\W £ - Wf\\ L 2+ 

(4.11) 

+\m - pIWv + M- = o(\\we\\v). 

We will use without further comments the elementary observation 

(4.12) || J f(x - s,k)g(s)ds\\ L 2( R 2 n) ^\\g\\ L i {Rn) \\f\\ L 2( R 2n) 

Lemma 4.1 (W ps W x ). V* G [0,T] 

Wm - W%t)\\ L 2 = 0( e-r^WWfWv ). 

Proof: Denote 

(4.13) h £ = W £ - Wf. 
Obviously, 

d t h + 2nk -d x h + T v \ = —T v ~ ^Vf, 



(4.14) 



h(x,k,0) = 0. 
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Since the Wigner equation has a bounded L 2 propagator, it suffices to 
bound in L 2 the rhs. Indeed, we have 
(4.15) 



T v - v Wf\\ L 2 = C\\JV(S){1 - e-l els2 )W{{X - S,K) 



sin(ireSK) 



dS\\ L 2 «C 



^0(e^)\\V(S)(l-e-^ s2 )\\ L i\m\\ L ^ 



f \V(S)\\S\ 2 dS + f \V(S)\dS 



\S\<e-i 



\S\^e-i 



0(s^)\\Wf 



li 2 



e^ + e^ f p~ e dp+ f p- 2 - e dp 



= 0{e~ l )\\W [ 
Asking that the remainder is small gives 

7 > 1, 



£ 7 + £ 7[-!(l-0)+l] +£ 7^ 



(4.16) 



1-6 



+ 1] > 1 <S> 7 > jf 



7^>1 * 7>tt^- 



So finally with the calibration 7 > the proof of Lemma 4.1 is complete 



Moreover: 
Lemma 4.2 (Wi « Vi G [0,T] 

(4.17) H^f(t)-pI(t)|| L2 = O(^-' 5 ||Ty e || i2 ). 

Proof: It is straightforward to check that 



(4.18) 



d t Wf - 2ttX ■ d K W{ + 2 / v (S)Wf(X - S, K) sm ^f- K ) dS = 0, 
dt% - 2irX ■ dxffi + 2vr / V(S)p%(X - S, K)S ■ KdS = 0, 



and therefore, if / = W £ \ — /5| 
(4.19) 



d t f - 2irX -d K f + 2f V(S)f(X - S, K) sin{n f J<) dS 
= 2 J V{SW 2 {X - S, K) (l - sin !:f;f ] ) ttS ■ KdS, 



■keS-K 

f(x,k,0) = W§-[W§\. 
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Using the Duhamel formula, 
(4.20) 

\\f(t)\\L^\m-m]\\ L 2+ 



+2T sup || J V(Sm(X - S,K,t) (l - Mn }:s S K K) ) *S • KdS\\ L2 . 

te[o,T] v 7 

Recall that the first term above is (relatively) small by assumption; we will 
work out the other term: 

To bound that, first of all observe that 



(4.21) 



sin(neS ■ K) 



neS-K 



1 



< Cmin{l,\ireS-K\} 



and therefore for any b < 
(4.22) 



/ v (S)MX - S, K) (l - s -^gjP) irS ■ KdS\\ L2 < 



<Cfe|| J \V(S)\-\S\^(X-S,K)\.\K\*dS\\ L2 + 

0<|S|<1 



+e 



1+261 



/ V{S)fF 2 (X-S : K)\K\HS\\ L 2 + 



K\S\<e b 



| / V(Sm(X-S,K)\SK\dS\\ L 2^ 

\S\>e» 



<C [ e^ b \\V(S)min{\S\\l}\\ L A\p £ 2 \\ H >+ I I V(S)\ \S\dSM\\ m < 

V \S\>e b J 



^ C ^e 1 +* b \\V 1 gr\\ L iM\\ B 2+ f^-'dpMWv J = 

= C(e 1 + 2b ||^|| ff2 +e- flb ||p||| H i ^C^+^ + e-^llWJHxa ) 

Calibrating the parameters is easy; we are given 6 £ (0, 1), 5 as in the 
statement of the Theorem and we need to find b 6 (—^,0) so that 



(4.23) 



1 + 26- <5 > 0, 

5-1 



5 > 



It should now be clear that the constraint on 5 in statement of the Theorem 
comes from the self-consistency check —\< < — | < 0. The proof of 



Lemma 4.2 is complete. 
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Now the last step is to show that p| « p\ ; but this follows by construction, 
since they satisfy the same equation (which has an L 2 -continuous propaga- 
tor), and | \p%{t) - p\(t)\ \ L 2 = \\W§ - [W^]\\ L 2. 



The proof of Theorem 2.2 is complete. 
4.3. Proof of Theorem 2.3\ One observes of course that, by construction, 

\E%0)\^M(t)-w%t)\\ L2 M\ L i. 

This will be useful in situations when we use a test-function from A n L 2 . 

Since in general there are cj) G A \ L 2 , we will also use a variation of the 
space of test-functions A introduced earlier: let 

(4.24) B = {<j> G C (R 2n ) | Sau P \^ K [<l>(x,k)]\dK + \\<l>\\ L 2<oo} 
equipped with the norm 

(4.25) IMb = IMU + IMb- 

In particular, it is clear by construction that B C A, A' C 5' , and ||/||l2 ^ 



Ib' ^ II/IIl 2 - Moreover, one easily sees the following 

Claim: There exists R > Osuchthat | \W £ {t)-pl (t)\ \ A > < R for all t G [0,T]. 
Proof of the claim: Observe that 

\{IA)a',a\ = I / Fk^K\f(%,k)}Tk^K[<t>{ x , k )]d x dK\ < 

< / sup\T k -> K [<i>(x,k)]\ f \T k ^ K [f(x,k)]\dxdK < 

< J s«p iJ^k^icl^x, k)]\dK sup J \T k -+K[f(x,k)]\dx, 

K x K x 



i.e. 



A , < sup/J : - fc ^[/(x,A;)](ix = snp/|/e- 2 ^^/(x,fc) ( iA:|dx^ ||/|| 



z 1 - 



To conclude we use the singular value decomposition of the initial data: 
if D £ = A m |n m )(u m |, then 

m 

W% = Y, A ™ [ e- 2mky ^ m {x + ^u^x - ^)dy, 

and 

W%t) =J2 X mf e-^uKx + f ~ £ {,t)dy, 

where of course each of u m (t), \\u m \\ L 2 = 1, satisfies the Schrodinger equa- 
tion (1.2). The point is that the singular values are time-independent and, 
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especially for positive operators, are controlled by the trace: 

D £ >0^\ m >0Vm => tr(D§) = Am = ||A m || J i. 

m 

We readily observe that 
sup J \Tk-*K[W e (x,k)]\dx = sup JY.\^mU £ m (x + ^f)u e m (x - ^f)\dx ^ tr{D e ), 

K x K x m 

and of course 

iipf(*)ib = iraui- 

The proof of the claim is complete; it is obvious that R ^ ^(-D^ + IH^ifllli 1 - 

We will also use the following: 
Claim: Let % G C°°(IR 2ri , [0,1]) be a cutoff function, X (z) = 1 V|*| < 1, 
x(z) = V|z| > 2, and denote <Pm(x, k) = X\m)- Then, 

n 

\\M\v* < rffS) Mn II^IU> 

and 

Urn \\4> M - 0IU = 0. 

M->-oo 

Proof of the claim: First of all, x(x,fc) G 5(M 2 ™) C A(R 2n ). Moreover, 
keeping in mind that ||0||l°° ^ H0lUi 



WIIl 2 ^ \m\L™ 



\ 



dxdk < M r 
y/T(n + l) 

|(x,fc)|<2M 



A- 



The second part, i.e. 



Jim ||0 M ~ 0| U = 

M->oo 



follows e.g. from the density of 5(M 2n ) in A jT3|. Additional information 
on the decay of is needed in order to estimate more precisely the rate of 
convergence. 

The proof of the claim is complete. 

It is clear now how, given a G A, we can select a family {0a/} C B such 
that ||0m — 0||.4^O as M^oo. Moreover set for brevity 

F{e) = C(e K \\W^\\ L , + \\[Wl]-W^\\ L ,) 

for some large enough C > 0. Now, recalling our remarks in the beginning 
of the section and the first claim, it follows that 

\{W £ - P \A)\ s: |(w e - P f,0 M )| + |(^ £ -pf,0M-0)| < 
(4.26) < F(e)||0 M || i2 + (1 + ||TV o £ || L i)||0 A/ - 0|U < 

< F ( £ )^gT) M " + C 1 + IMIia)||0m - 0IU Ve,M. 
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The proof is completed by selecting, for each value of e a "small enough" 
M, so that M = M(e) -> oo, and 

limF(s)\\4> M ( e) \\ L 2 = 0. 

e-»0 

In the statement of the Theorem, the calibration M = F(e)^^i is used; in 

n 



that case it follows that II^m(£)IIl 2 ^ 



The proof is complete. 



Vr(n+i) 



5. Examples 

This section is devoted in understanding better the non-concentration 
properties that are in the heart of our main results. There are two versions; 
we isolate them below as "Property 1", a condition to be checked for a 



given problem that is sufficient for Theorem 2.2 to apply to it, and the 



stronger "Property 2", which is sufficient for Theorem 2.3 and its Corollary 



2.4 Concrete examples of natural problems that satisfy either one or both 
of them are constructed; it must be noted that it's rather easy to satisfy 
Property 1 to generic initial data and repulsive singularities. Asking that 
Property 2 holds is quite restrictive, and it seems that specially prepared 
initial data are needed for that. 

We are only concerned here with examples that give rise to ill-posed 
classical problems. (Otherwise the existing theory for C 1 potentials has 
no problem to provide a complete description of the semiclassical limit). 



5.1. Examples for Theorem 2.2\ Here we are concerned with Cauchy 



problems for the Liouville equation, 

dtpl + 2irk ■ d x p\ - ^d x V -d k p\ = 0, 

pf(t = 0)=W§, 
which have the following 

Property 1. There exist T > 0, 5 G (0, ^) and [Wfi] G H 2 (R 2n ) such that 
(5-1) \\m]-WS\\ L 2=o(\\WS\\ L 2) 
and the solution of 

dtpl + 2n k • dxP £ 2 - ^d x V -dkpl = 0, 
^ = 0) = ^], 

satisfies 

(5.3) \\pl(t)\\ H 2 = 0(e- S \\W^\\ L 2), uniformly for t G [0,T]. 
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We treat k, 7 as given parameters; of course when checking the property in 



the context of Theorem 2.3 they are controlled by the statement of Theorem 



Lemma 5.1. Let 6 G (0, 1), 5 G (0, n = 2, V(x) = -\x\ 1+e with an 

appropriate smooth cutoff outside {\x\ < 2}. Assume moreover that there 
exists a function of compact support w G ff 2 (IR 4 ) n L°° , such that 

(5.4) ^, fc ) = i 3-M^p),^^), 

with (zijZz) = —L(z^,z^) (i.e. the wavepacket is "shot towards zero"). 
\\w\\ H 2, \\w\\l°°, \z\ have to be bounded uniformly in e, \\w\\ H 2, \\w\\l°°, \z\ = 
0(1); other than that w and z can be allowed to depend on e. 

The small parameters involved S x , 5 k , R = o(l) are calibrated as follows: 

1 

1 \ \ s 

R=\log\ 



as in equation (6.26), and 

8- 1 (R + 5 k ) = o(l), 

(5.5) 

{5 x + 5 k )~ 2 {5 k + R)- 2 = 0{e-2 



Then Property [7] is satisfied. 
Remarks: 

• In fact it follows from the Lemma that property [T] holds for any 
5 €(0,^5). 

• The difficulties arise from the non-smootheness at zero, and the de- 
tails of a smooth behaviour away from zero are irrelevant here. So we 
will work for a small enough time, before any trajectory starting in 
a neighbourhood of zero reaches the support of the cutoff function. 
Therefore, to keep the presentation simple, we will not introduce any 
explicit treatment of the cutoff function. This approach is followed 
in the sequel as well. 

• An example of a scaling satisfying the above constraints is 8 X = 
Vh = VR. 

Proof: Assume without loss of generality that z = (0, — L, 0,1). Consider 
<fi : K -» [0, 1] to be a C°° function such that 

, s <t>(x) = 0, \x\<l, 

1 ' <f)(x) = 1, \x\ > 1. 

The modified initial data [W§] is selected as 



(5.7) [W£](x,k) = Wt(x,k)cf> 



\xi\ 
R + 2L5 k 
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The key claim is that, for T small enough but independent from e, 

(J supppi(i) n {\x\ < R} = 
te[o,T] / 

Indeed [Wq] has two disjointly supported components. The calibration of 
the parameters and the cut-off is such that, if [Wq] was propagated by 
a Liouville equation corresponding to V(x) = 0, each of the components 
would stay on either side of the strip {|xi| < R} (see Figure [2]). The idea 
is that since there is a potential driving the flow away from x = 0, the 
claim follows. (In fact it is slightly more complicated for this case, since the 
potential depends on x% as well: one easily gets an O(l) upper bound for 
the time it takes for the trajectories to reach x% = 0, and then checks that 
in this longer time, x\ doesn't have the time to change sign. The difference 
in an 0(1) factor which we absorb in C below). 




Figure 2. Any trajectory leaving the support of [Wq] has 
an initial velocity in a small cone around k = (0,1). If we 
were in free space, there wouldn't be enough space for suf- 
ficient movement in x\ to reach {|xi| < R}. Therefore the 
solution is the sum of two components, supported on either 
side of {\x\\ < R} = 0. The presence of a repulsive potential 
in many cases can be easily factored in this construction. 



~7 ~7 / II \ 

Therefore, in this problem we can take V e ff(x) = V (x)<j> ( R+cL8 k ) • Now 
for condition (5.3), making use of observation 6.6 (enlarging the cut-off area 



in fact strengthens observation 6.6 i.e. substituting R + 2L5k for R is 
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painless) and Lemma 6.4 it follows that 



(5.8) 



\pm\\H* < \\m]\\H^-^c\\ws\\ H 2e-H5 k + Ry 

= C\\W$\\ L 2(5 x + 5 k )- 2 (5 k + R)- 2 e-l, 



and therefore - using the scaling of equation (5.5) - | |p|(*)ll.ff 2 = 0( £ 5 | I I II 2 )- 

Now for n - it;,;: 



^ 5; 



-2ns-2n 



J 



W 



' (xi,x2) (fci— l,fea) -\|2 



\xi\<C(R+6 k ) 



S x ' 



)\ 2 dxdk 



(5.9) 



f \w(xi, X2, ki — 1, k2)\ 2 dxdk ^ 



^C7||^||| 2 ||u;||| oo( 5- 1 (fi + 4) 
The same idea can be applied to different configurations: 



Lemma 5.2. Consider the setup of Lemma 5.1 with z = (0,-2,0,1) and 
the only difference that 

(5.10) V(x) = -\x 1 \ l+6 i>(x 1 )ip(x 2 ). 

where ip is a smooth cutoff function, ip = 1 — 4> £ S(M) ((p was defined in 
equation (5.6)). 

Then Property^ holds. 

Proof: The geometry is essentially the same as before, and it is clear that 
the obvious adaptation of observation |6.6[ i.e. the one with 

holds. 



5.2. Examples for Theorem 2.3 In relation to Theorem 2.3, we intro- 
duce 

Property 2. Given V(x), Wq with \ \Wq\\ L 2 = o(e~ K ), there exist T > 0, 
5 G (0, 2^g) and M e H 2 {R 2n ) such that 

(5.H) \\m]-w^\\ L 2 = o(i) 

and the solution of 

dtpl + 2irk ■ d x p\ - ^d x V -d k p% = 0, 
pl(t = Q) = [W^ 



(5.12) 
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x2 



21! 



x1 



Figure 3. The construction of Lemma 5.1 holds, with ob 



vious adjustments, to the example of Lemma 5.2 



satisfies 
(5.13) 



0(e 



—5 1 



uniformly for t G [0, T]. 



As was already seen in section 



4.1 



without loss of generality we can check Property [2] in the context of Theorem 




-Z7T- 



\ L 2 = o(l). So 



2.1 considering initial data of the form /g , and not the smoothed version that 



appears there. (More generally, working with compactly supported initial 
data - and hence removing the smoothing - might be very helpful whenever 
we check Properties [T] [2] In any case, we expect that the smoothing can 
be removed without loss of generality, just like in this case, whenever the 
initial data decays rapidly and ||Wg||#i = o(e~2)). 



Lemma 5.3. Let n = 1, 6 € (0,1), V{x) = — \x\ 1+e with an appropriate 
smooth cutoff outside {\x\ < 2}. Assume moreover that w(x, k) E H 2 nL°°n 
L 1 , suppw C {|x|, \k\ < 1}, and 



(5.14) 

Finally, set R = 
Then Property 



A: 



1 w(— — ) 



n is satisfied if 5k 



a = C"Rs- 
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Proof: We will cutoff a strip of the form {|x| ^ CR}, and show that in fact 
this suffices. 

As is illustrated in Figure 5.2 the preimage under the flow of {\x\ < 

R A \k\ < 5k} is contained between the level sets {2ir 2 k 2 — V (x) = V (0)} 

and {2n 2 k 2 + V(x) = V(R)}. 

A big help with the algebra will be to ap prox imate V (x) = — \x\ 1+6 + 
0(ez). (7 is as in the statement of Theorem 2.2). Since all the other small 
parameters are (negative) powers of log(^), O(e^) is negligible everywhere. 
(The justification is that without loss of generality we can localize the prob- 



lem on a compact set, and then V(x) £ W '°° uniformly in e, see Lemma 6.5 



The conclusion follows by a standard observation on mollifiers, see Lemma 



6.2). 



It is easily seen that the needed length is x*, defined by 

2ir 2 5 2 k + V (x*) = V{R) 

or 

\ x ^+o = 27r 2 5 2 + R 1+e + 0(ei). 

Here is where 5k = C'R^ comes from; making that scaling we get 
x ^ = ((1 + 2ir 2 )R 1+d + 0(e?))TTfl = C'R. 

So far we have ensured by construction that, if 5k = R~^~ , (there is an 
0(1) constant C so that) if 

(5-15) [Wl]{x,k) = f Q {x,k)<l ) {-^ i ), 



its propagation under equation (5.12) will never enter {|x| < R}. 
Now for the approximation error: 

u-m]\\h^c5-\ 2 j \ W (f x ,f- k )\ 2 dxdk = 

\x\<CR 

(5.16) 

= C6~% 1 / \w(x, k)\ 2 dxdk < C5~ 2 RN L . 

\A<cf x 

This gives the constraint R~ = o(5 x ) for the error to be small, which is 

1-9 

satisfied e.g. by our earlier selection 5 X = CR 5 . 

Again, since all the small parameters are powers of log(^), ||[Wq]||#2 = 
o(e~5) follows automatically for any 5 £ (0, ^g) ■ 

The proof is complete. 

The following is a somewhat artificial example, but it highlights an inter- 
esting behaviour. 
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Figure 4. Schematic for Lemma 5.3 cutting off a strip in 
x, {\x\ < x*}, suffices to cut off the preimage of {|x| < R} 



Lemma 5.4. Assume n = 2, V as in equation (5.10). Moreover, there is a 



function f G 5(R 4 ) ; f f(x, k)dxdk = 1, such that 
(5-17) 

Denote by f 
(5.18) 
Now take 

w§(x, k) = ^ k ( cl / r- 2 £" +2 , ^f 1 



supp f C {|x| < 1} x {|&| < 1 A ki > 0}. 
f(xi,X2,ki,k 2 ) = f(-x 1 ,x 2 ,-k 1 ,k 2 ) 



+ 



(5.19) 



+C2 ->( — s» — 'Si—- 



If $k = C'&z = C" y R, this problem has Property^ 

Proof: Two things should be obvious by construction: that W§ 5(x + 
(0, 2), k — (0, 1)), and that the propagation of Wg under equation ( |2.8[ ), i.e. 
p\ , is never supported inside < i?}. In other words, for this specifically 
constructed data, Wq = [Wq]. 

The scaling of 6 X , 5^ in this case controls simply the rate of concentration, 
and therefore ||Wq||#2. The result follows. 
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6. Auxiliary results 

Lemma 6.1 {2 nd order derivatives equations for the Liouville equation). 

Consider the Cauchy problem for the Liouville equation with potential V(x) 
on R n , 

d t f + 2nk • d x f - ±d x V(x) ■ d k f = 0, 

(6.1) 

f(t = 0) = /„. 

There are constants C\,Ci > depending only on n such that 

tC 2 sup \\d°V(x)\\ L oc 

\\f(t)\\ H ^c ie ^ \\M\ H 2. 

Proof: The proof follows readily with the method of characteristics. 
It suffices to observe that if 

(6.2) ±i{t) = 2irki(t), k{t) = -^d Xi V(x(t)), 

z(t) = f(x(t),k(t)), 
z Xi {t) = d x J(x(t),k(t)),z XiXj (t) = d XiXj f(x(t),k(t)), 

kj (t) = dkikj 
z Xi kj{t) = d Xikj f{x(t),k(t)), 

it follows that 



(6 ' 3) z ki (t) = d k J(x(t),k(t)),z kikj (t) = d kik J(x(t),k(t)), 



(6.4) z = 

(6-5) z Xi (t) = ^Y J d *i*mV{x(t)) z Xikm (t), 

m 

(6.6) z kt (t) = -27Tz Xi (t), 

( 6 - 7 )= «- [ 9 xjx m Vz Xikm (t) + d XiXm Vz Xjkm (t) + d XiXjXm Vz km (t)] 



2vr 

m 



(6.8) z ktkj {t) = -2ir(z Xjki + z Xikj ), 

(6.9) z kiXj {t) = -2nz XiXj (t) + ^^d XmXj V{x{t))z krnki {t). 

m 

The result now follows by the Gronwall inequality. 
Lemma 6.2. Consider a function f e W 1 ' 00 , i.e. sup \ \d x f\\L°° < °o- 

|o|<l 

Then, if f ' = f - J f e i f(x l )dx* , we have 

(6.10) ||/ -711^=0(11/11^00^)- 
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Proof: Take any £ G (0, \). Now we have 

!/(•>•) -/<•>•>!- I 2 ]2 ' 2t 



27 



IV f e - Z7r ^[f(x-x')-f(x)}dx> 



2 \ 2 

'/ 



/ e 2w v [f(x — x)' — f{x)]dx' 



|x'|<r?2- 



+ 0(77°°) «S 



3.111 



2 \ 2 



J e- 27r ^t\x'\dx' + 0(11/11^°°) 



|x'|<??2 C 



= 0(^2^ 
The proof is complete. 

Remark: A sharper version is used in subsection 4.1, namely 
/ \ 



w ^ f e-**\vny\dy + 0(\\f\\L°°V c 
\y\<v~ c 

= 0(||/|| W l,oo^?). 



\f(x)-f(x)\^C 



sup \d£f\ 

\ \x - x'\ < ri2~^ 



2 \ 2 
1j 



f e 2n v \x'\dx'+ 



|x'|<r?2" 



+0(\\f\\L-T, C 



The following observations are used in Section [5] 

Observation 6.3 (Locality of the Liouville equation). For a Liouville equa- 
tion with initial data of compact support pq and a C 1,1 potential, interchang- 
ing the potential with any one that coincides with it on the "path" of the 
solution, 



3.12) 



§ = (J 4>t(supppo) 
te[o,T] 



does not change the solution. 
Also, 

Lemma 6.4 (An H 2 estimate). Consider equation (2.8) with initial data 



9 

Pq = [Wq] S H . Assume that there is a function V e ff such that 
(6.13) 



~7 ~7 

V = Veff on S, 



2N 
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where § is defined as in (6.12\ ), and 

(6.14) sup \\d£V e ff{x)\\ L ™ =o(log(e- 

|A|<3 V V 

Then, 

(6.15) HpIH^ =0(e-l\\W^\\ H 2 ). 

~7 ~7 

Proof: The proof consists of using V e ff in place of V , making use of ob- 



servation 6.3 and then applying directly Theorem 6.1 



We will also use the following 

Lemma 6.5. It is easy to observe that, ifW G L°° and suppW^ is compact, 
it follows that 

(6.16) 3M > : Vt G [0,T], e > : |J supp p\(t) C {|(x, fe)| < M}. 

te[o,r] 



Observation 6.6 (V e // and R). Let V(x) 
version, according to equation (2.6). 



x | 1+9 . V is a mollified 



Consider (f> : R n -> [0, 1] to be a C°° function such that 

x) = 0, |x| < 2> 



3.17) 



<t>(x) 



1, a; > 1, 



and ^10. is clear that such a function exists; this is of course an 

arbitrary requirement, but one that allows us not to carry the cutoff function 
(j) to other results. 

Then, setting R ^ {log (e -1 )) 5 , R = o(l), it follows that 

(6.18) 
satisfies 

(6.19) 

while, of course, 

(6.20) V eff = V on R n \ {\x\ < R}. 



Veff(x) -=<Kjj) V(X), 



sup \\d£Ve ff (x)\\ L ~ =0(log(e-l) ), 

|A|<3 V ' 



Proof: Observe that 1 

observations 
generality, 



6.5 



and 



6.3 



(j>(x) = ip{x) G 5(M n ). Moreover, making use of 

~7 

we can restrict V to {\x\ < M} without loss of 



3.21) 



v\x)^ i7(*)(i-0(JL)). 
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Denote 

K // (x)=y(x)(l-0(^))0(l) 



V(x) (V(dr)-^(S)) 



(6.22) 
Since 

(6-23) \\d£f\\ L oo < \\d£f\\ L oo, 

it suffices to work with 

(6-24) (*(£)-*(*)) 

Set x' = J^; then <9 X . = R~ 1 d x i. Now we have 
(6.25) 

ll#M 14 * ^(^r)-^(5))l|£- 
= J R e - 2 ||^|x| 1 + e (^(^x)-^)) 



L 3 



< R 2 SUP ||3^X 1+9 || LOC( rl 3M]) SUP \\d£ (l/j [jtfX) - i>(xj) ||l/x>(]R™) < 

|A|<3 2 ' R |A|<3 

< {2M) l + e U\\ w ,, ao R- ? \ 

Recall that the constant M is chosen so that any trajectory leaving the 
support of the initial data doesn't not escape {|x| < A/} for t £ [0, T}. That 
is, to properly quantify it, one needs to consider initial data of compact sup- 
port associated with the Liouville equation, and a time-scale T. Assuming 
that the initial data is supported in {\(x, k)\ < Ro}, it is easy to ch eck that 
one can set M = R + 1 + RqT + | |W(x)| \ L °°T 2 . So now for ( |6.19| ) to hold 
it suffices that 



2 1+e (1 + Rq + R T + \\VV(x)\\ L -I a ) 1+t 'io d _ 3 
C 



R~ 3 ^ log(e 



-In 



2 

The constant C comes from the O(-) of equation ( 6.19| ); it is clear that it can 
be chosen so that the constraint finally becomes R ^ (log(s~ 1 ))^s , which is 
satisfied by choosing 

(6.26) R=(log(e- 1 ))-l 

The following is contained in [T3]; we include here a brief mention for 
completeness: 

Lemma 6.7 (Density matrices). Let fj, be a probability measure on M. 2n . 
Then 



2 x n 



W$(x, k) = 1^1 J e~ 27r - * dn(x', k') 

is the Wigner function of a density matrix, i.e. the corresponding operator 
D £ is a positive trace-class operator with tr(D £ ) = 1. 
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Proof: We know that tr(D e ) = J W £ dxdk = 1. Let us now look at positiv- 
ity. To that end, observe that the integral kernel 

2y J e 2 W i^(*-») e -f[( a! -xo) a +(v-*o) 3 ] d/1 ( a . 0jA!o ) 

x ,k 

is the kernel of a positive operator. Indeed: 
fK £ (x,y)u(x)u(y)dxdy = / \( e - 2 ^ x+ T n)| 2 d/i(x , k ) > 0. 

x ,k 

The proof is complete by observing that the Wigner function correspond- 
ing to the kernel p £ is 

J e -2mky K e {x+ ^^_^_ )dy = w e_ 
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